Draft version February 16, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



O 

o 

D 



6 



> 

o 

(N 

(N 
O 

a\ 
o 



SPECTRAL STATES AND EVOLUTION OF ULTRALUMINOUS X-RAY SOURCES 

HuA Feng^'^ and Philip Kaaret^ 

Draft version February 16, 2009 

ABSTRACT 

We examined spectral evolution in ultraluminous X-ray sources (ULXs) with apparent luminosities 
of about lO""^ ergs s~^. Based on new results in this paper and those reported in the literature, two 
common spectral behaviors were found. Some ULXs in starburst galaxies have varying luminosity (L) 
but remain in the hard state with power-law spectra and a constant, hard photon index (F). Other 
ULXs, such as NGC 5204 X-1, show a correlation between L and F. We interpret this L — F correlated 
phase as an intermediate state with hybrid properties from the thermal dominant and steep power-law 
states. When the spectra of NGC 5204 X-1 are fitted with a multicolor disk blackbody plus power- 
law model, the X-ray luminosity increases with the effective temperature of the accretion disk in a 
manner similar to that found in stellar-mass black hole X-ray binaries, suggesting that the emission 
arises from an accretion disk. The luminosity, disk size, and temperature suggest that NGC 5204 X-1 
harbors a compact object more massive than stellar-mass black holes. In contrast, the disk model in 
IC 342 X-1 is ruled out because the luminosity decreases as the temperature increases; sources with 
such behaviors may represent a class of objects with super-Eddington accretion. Also, we report a 
peculiar soft spectral feature from IC 342 X-2 and variability on a time scale of 20 ks from Holmberg 
II X-1. More observations are needed to test these results. 

Subject headings: black hole physics — accretion, accretion disks - X-rays: binaries - X-rays: indi- 
vidual (NGC 5204 X-1, Holmberg II X-1, IC 342 X-1 and X-2, the Antennae X-11, 
X-16, X-42, and X-44, NGC 1313 X-1 and X-2, M82 X41.4+60 and X42.3+59) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are point-like, 
nonnuclear X-ray objects with apparent luminosities (in- 
ferred assuming isotropic emission) above the Edding- 
ton limit of a stellar-mass black hole, typically 3 x 
gj-gs g-i fQj. a, 20 Mq black hole. Those with 
rapid variability must be a.ccreting compact objects (e.g . , 
IStrohmaver fc Mushot^l2003l : [St"rohmaver et aLll2007h . 



Their high luminosities suggest they might be powered 
by accretion onto black holes with masses in the range of 
10^-10'* Mq, i.e. the missin g population of inte rmediate- 
mass black holes (IMBHs: iColbert fc Mushotzkv 19991 : 
iMakishima et al.l 120001 : iKaaret et all l2001l ). Alterna- 
tively, they could be a special class of stellar-mass 
black hole binaries, different from most ones in our 



sion (|Fabrika & Meschervakovl 


2001: King et all 120011: 


Watarai et al. 120011; iBegelman 


120021 IBegelman et al.l 


20061; iPoutanen et al.ll2007D. Although their uhvsical na- 



ture is still unclear, ULXs have shown diver se phenomena 
and may be a hete rogeneous population (jSwartz et al.l 
[200l iFeng fc Kaaret 2005. 2006b). 

An important argument that supports the IMBH sce- 
nario is the discovery of soft excesses in the energy spec- 
tra of ULXs, whic h can be interpreted as emission from 
an accretion disk (i Kaaret et al.l [20031 ). At a fixed frac- 
tion of Eddington luminosity, the disk luminosity is pro- 
portional to the black hole mass, while its inner tem- 
perature is scaled with the mass to a power of — 4 
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(|Makishima et al.ll2000l) . Disks in U LXs have inner tcm- 
peratures in the range of 0.1-0.4 keV ( Miller et al. 2004^; 
iFeng fc Kaaret l2005t iStobbart et al.( 12006 V cooler but 
more luminous than in stellar-mass black holes, sug- 
gesting that the emission arises from accretion onto 
IMBHs. However, th e cool accretion disk mo del has 
aroused broad dispute. iGoncalves fc Sorial (|2OO60 argued 
that the soft excess could be a soft deficit depend ing on 
how th e power-law continuum is modeled. Berghea et al.l 
(|2008f ) found similar soft excesses exist in both ULXs 
and less luminous obj ects, i i nplyin g tha t they are not a 
signat ure of IMBHs. iKind (|2004D and IPoutanen et all 
(|2007^ suggested that the cool thermal emission is as- 
sociated with massive outflow in super-Eddington ac- 
cretion rather than from the inner disk. Due to the 
quality of the data, some spectra could be equally well 
fitted by other models that do not need the presence 
of IMBHs ( Stobbar t et al.ll2006l;lvierdavanti et al.ll200l 
iMizuno eF al. 2007). Study of multiple spectra obtained 
at different luminosities from individual ULXs enables 
one to test if the spectral parameters evolve as expected 
in the models. This provides a means to resolve the ambi- 
guities in fitting individual spectra and determine which 
models correctly describe the physical behavior of ULXs. 

Due to the nature of blackbody emission, an opti- 
cally thick accretion disk should follow the L oc T* 
relation, where L is luminosity and T is the tempera- 
ture. This has been widely v er ified in black ho l e bina- 
rics (Gicrlihsk i fc Donell200l . iFeng fc Kaaretl (|2007aD 
found that the thermal temperature of NGC 1313 X-2 
failed to follow the L oc relation while the luminos- 
ity changed by a factor of 10, providing solid evidence 
against the disk model. However, it is interesting that 
the spectral evolution of NGC 1313 X-2 shows a pattern 
totally different from what was observed in NGC 1313 
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X-1 (jFeng fc Kaarelll2006aD . The cool disk model, al- 
beit not favored statistically and proved wrong in one 
case, still could be correct in some other ULXs. Finding 
spectral evolution with a pattern of L cx T'^ would be 
evidence in favor of the cool disk model. 

Spectral evolution could also help us understand the 
spectral states of ULXs. Black hole binaries may show 
similar spectral features in different emission states. 
Thus, without timing or multiwavelength information, 
the spectral state often cannot be unambiguously deter- 
mined with a snapshot observation. Multiple observa- 
tions showing an evolution pattern would be more con- 
vincing in determining the spectral state. 

We investigated spectral evolution in eight ULXs, NGC 
5204 X-1, Holmberg II X-1, IC 342 X-1 and X-2, and An- 
tennae X-11, X-16, X-42, and X-44, using archival and 
proprietary XMM-Newton observations, as well as Chan- 
dra data if necessary. These sources have peak luminosi- 
ties over 10"*" ergs and strong long-term variability, 
thus are good candidates for spectral evolution search. 
Some ULXs like NGC 5408 X-1, though of great inter- 
est, have shown roughly constant spectral parameters in 
the literature and are not included in the paper. The 
details of observations and data analysis are described in 
§ m In § [3l we discuss our results together with results 
on the spectral evolution of other ULXs with reported in 
the literature. A brief conclusion of the paper is in § 

2. OBSERVATIONS AND DATA ANALYSIS 

We used SAS 7.1.0 with calibration files current as of 
2008 July for X MM datSL reduction. New event files were 
created from Observation Data Format (ODF) files for 
each exposure. Multiple exposures made with the same 
Charge Coupled Device (CCD) in the same observation 
were combined into a single spectrum (data from dif- 
ferent CCDs are not combined). Time intervals were 
excluded if flares occurred in the 10-15 keV lightcurve 
of the whole CCD. Energy spectra were extracted from 
events with FLAG equal and PATTERN no more than 
4 for FN or 12 for MOS. For each observation, we used 
the same set of source and background regions in sky co- 
ordinates to extract spectra from different CCDs. The 
source extraction region was defined as a circle centered 
on the source with a radius of 32"; if there were bad 
pixels or CCD gaps within the region, we shrunk the cir- 
cle until they were excluded. Therefore, the extraction 
regions could vary for different observations. We note 
that the size of the source extraction region does not 
affect the inferred flux. The background regions were 
selected to be near the source, off the readout column, 
and on the same CCD chip as the source. Energy chan- 
nels of each spectrum were grouped to give a minimum 
of 25 counts per bin. Spectral fits were performed us- 
ing XSPEC 12.4 available in HEASOFT 6.4. The errors 
on the flux and luminosity (absorbed and unabsorbed 
fluxes) were calculated in XSPEC 12.5 with the cf lux 
model. Background subtracted spectra from MOS and 
PN were fitted simultaneously with the same model in 
the energy range 0.2-10 keV. Constants to account for 
possible discrepancy between MOS and PN normaliza- 
tions were not necessary, thanks to the identical extrac- 
tion regions and the relative dimness of sources. The 
absorption column density in our Galaxy along the line 
of sight to each source was obtained from the LAB map 



TABLE 1 

XMM OBSERVATIONS OF NGC 5204 X-1 



Good Exposures (ks) 



No. 


ObsID 


Date 


PN 


MOSl 


MUS2 


1 


0142770101 


2003-01-06 


15.3 


18.5 


18.5 


2 


0142770301 


2003-04-25 


3.8 


7.6 


7.6 


3 


0150650301 


2003-05-01 


4.9 


7.9 


7.9 


4 


0405690101 


2006-11-15 


9.5 


17.4 


17.6 


5 


0405690201 


2006-11-19 


29.0 


39.7 


40.7 


6 


0405690501 


2006-11-25 


20.4 


29.9 


30.0 


Note. 


— Good Exposures are 


effective 


exposures after 



background screening. 



(jKalberla et al.ll2005D using the tool on HEASARC, 
and was set as a lower bound of the total column density 
in the absorption model. Other parameters specific to 
each observation are described individually below. All 
quoted errors in the tables are of the la level in order to 
be consistent with plots. Since the errors on the lumi- 
nosity are usually smaller than on the temperature and 
photon index, we treat the luminosity as an independent 
variable and temperature or photon index as a depen- 
dent variable in the curve fitting, and take the errors on 
the dependent variable only to calculate x^- 

2.1. NGC 5204 X-1 

NGC 5204 X-1 has recently be en observed multiple 
times with XMM and Chandra ()Roberts et all l2005l 
l2006l) . A monitoring program with ten Chandra snap- 
shots revealed a flux change by a fa ctor of 5 at timescales 
of a few days (|Roberts et al.ll2006l ). By aligning Chan- 
dra and HST images, an optical counterpart to the X- 
ray source was i dentified, which appears like an early 
type supergiant (|Liu et al.ll2004l ). Here, we re-analyzed 
all six archival XMM observations of the source listed 
in Table [T] with observation IDs, dates, and "good expo- 
sures" that indicate effective exposure times after remov- 
ing background flares. The first t wo observations have 
been reported in the literatu re (e.g., [R oberts et al."2005|; 
. Feng fc Kaa ret 2005; Stob bart et alT .2006: Winter et a] 
l2006t fVierdavanti et al.ll2006i r while the later four have 
never been published. Our source extraction region has 
a radius of 27.5" for the first two observations, 20" 
for the 3rd, and 25" for the last three, respectively. 
We adopted a dis tance of 4.3 Mpc to the host galaxy 
(jTuUv et al.lll99^ . and a Galactic absorption column 
density of 0.174 x 10^^ cm"^. 

We tried different models to flt the data. All spec- 
tral parameters, observed fluxes /x and unabsorbed lu- 
minosities Lx in the 0.3-10 keV band, x^, and degrees 
of freedoms are listed in Table [2] If all model compo- 
nents are physical (with a finite integration), a bolomet- 
ric luminosity Lboi is also calculated by integrating the 
model (without absorption) in the energy range of 0.01- 
100 keV. First, we tried an absorbed power-law model 
(powerlaw in XSPEC), as the first-order approximation 
of the spectral shape. A cool disk emission component 
(diskbb in XSPEC) was then added, which improved 
the fits significantly. To replace the power-law compo- 
nent with a physical model, we tried a Comptonization 
model (comptt in XSPEC), with and without a disk com- 
ponent. When the disk component is added, the disk in- 
ner temperature Tin and the Comptonization seed pho- 
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TABLE 2 

Best-fit parameters of NGC 5204 X-1 



No. 


Nil 


r/r 


Npi^/Nc 


Tc Tin /To RinV^ 


/x 


Lx 


-C-bol X^/dof 


(1) 


(2) 


(3) 


(4) 


(5) (6) (7) 


(8) 


(9) 


(10) (11) 


Model: wabs*powerlaw 



1 


n 5Q+0 03 

'^■"'^-0.03 


2.07 


2 


'^•yj^-0.06 


2.29 


3 


1 ,0+0.06 
'^■'^■^-0.06 


2.42 


4 


1 41+0.04 
^•^^-0.04 


2.66 


5 


1 40+0-03 
^■^'^-0.03 


2.561 


6 


n oq+0.03 


2.251 



+0.02 
-0.02 
+0.04 
-0.04 
+0.04 
-0.04 
+0.02 
-0.02 
+0.016 
-0.016 
+0.020 
-0.018 



3.58 
5.60 
7.06 
10.06 
8.35 
4.95 



+0.06 
0.06 
+0.17 
0.17 
+0.20 
0.19 
+0.18 
0.18 
+0.10 
0.10 
+0.07 
0.07 



1.62 



+0.02 



-0.02 
^+0.04 
-0.04 
1+0.04 
-0.04 
1+0.03 
*-0.03 
9 9QS+0.015 

1 8fifi+0018 
l-8o6_o.oi8 



2.02; 
2.241 
2.1 



n 4on + 0.005 

3+0.015 
-0.014 



0.619^ 



0.761 
1.08 
0.892 
0.553 



+0.019 
0.017 
+0.02 
0.02 
+0.011 
0.011 
+0.006 
0.006 



513.1/461 
273.8/219 
288.3/262 
593.5/566 
926.8/817 
807.0/658 



Model 



wabs(powerlaw + diskbb) 



1 


56+'''05 


1 Q1+0.04 

-^•y-^-0.05 


9QO + 0.16 
^•^'-'-0.19 


26+0-03 


75+0-16 
'0-0.15 


1 65+0-02 
l-°°-0.02 


430+0-008 

U.'iOU_g ggg 


491.8/459 


2 


71+0-14 

'^•'1-0.10 


1 09+010 
^■y^-0.10 


q c+0.5 


29+0-03 
"•^^-0.03 


1 01+0-33 


9 1 1 +0.04 
^■-^-^-0.04 


50+0.03 

'-'■0^-0.02 


249.1/217 


3 


n 09+0.15 

0-»^_0.09 


2 10+0-12 

^-I'-'-O.lO 


4 6+0-9 


qq+0.03 

U.OO_g gg 


72+0-24 
'^-0.12 


2 -^0+0-04 
^-■'"-0.04 


PI fio+0.03 


277.1/260 


4 


1 06+0-12 
^•^°-0.07 


9 .9+0. 08 
^-^^-0.06 


7 0+0-9 

'■"-0.6 


36+0-02 
"■•^"-0.03 


64+0-12 




90+0 04 

'-'■y^-0.04 


578.3/564 


5 


94+0-07 

'^•y^-0.04 


9 ,0+0.05 
^■18-0.04 


A S+0-4 
4-8-0.3 


PI qpr9+0.009 


p. -0+0. 06 
'J- '■J-0.04 


2.348lg:g- 


p, 799 + 0.017 

U.(ZZ_ggig 


837.0/815 


6 


63+0-04 
'^•°'^-0.06 


1 ot:+0.04 
l-^''-0.06 


2 96+0-18 


907+0.018 

U.Za( _g g-^g 


90+0-0'7 

'j-y'^-o.io 


1.96«:«2 


n c:99+0.009 

U.OZZ_g ggg 


697.7/656 



Model: wabs*comptt 



1 


0.17<0-2i 


5 61+0-08 
^•°l-0.54 


4 0+0-3 


9 sS+0-19 
2-OO-0.23 


116+0-003 

U.liD_g ggg 


1 61+0-02 
1-O1-0.02 


074 + 0.005 
U.jr4_g gg4 


0.44 


492.7/459 


2 


0.17<0-25 


1^:1+0.012 

U.101_g ggg 


n 149+0012 

>J.14^_g.003 


[100] 


0,4, +0.003 

U-l*l-0.005 


2 04+0-04 
^■'J^-0.04 


n 471+0.022 
'J-^'1-0.008 


0.67 


243.2/218 


3 


29+0-14 
'^•^y-0.07 


,9r-+0.009 

U.iZO_g gg7 


f. ,79+0.024 

'-'■1 '^-0.016 


[100] 


142+0-005 
U.14^_ggg3 


9 90+0.04 

Z.ZO_g g4 


55+0-02 
'^■"^"'-0.02 


0.73 


275.9/261 


4 


26+0-08 


4 8+0-6 


9.+2- 


2 4+0-^7 


152+0-007 

U.10Z_g ggg 


9 p.9 + 0.03 
Z.DZ_g gg 


0.628tlfil 


0.67 


570.5/564 


5 


,03+0.014 

U.iyj_g gj4 


0.0983+r»JI 


160+0-005 

U.iDU_g gg- 


[100] 


1647+0-0014 

U.iD4r_gggi3 


9 q9c+0.015 
Z.OZD_g gj^g 


0.542i0.o;5l 


0.68 


828.7/816 


6 


0.17<0-2i 


0-167t°:»I 


1315+0-0065 

U.1010_g ggj^2 


[100] 


1359+0-0016 

U.iOOy_o ggg2 


1 SSQ+0-021 

i.»»y_ggig 


0.438l«:r3 


0.64 


728.4/657 



Model: wabs(comptt + diskbb) 



1 


29+0-05 

'j-^y-o.os 


10.3+12 


9 c: + 0.3 
^■^-0.3 


1 50+0.17 

i-3y-o.i5 


2 


o.5ilS:iJ 


20+0-05 
'^-^'-'-0.03 


0-08212:°?^ 


[100] 


3 


o.5ot«:i; 


PI 109+0.024 

U.iCiZ_g gjg 


0.12lH^ 


[100] 


4 


52+0-06 
"•''^-0.07 


5.2+?-? 




9 9 + 0.9 
^■^-0.4 


5 


56+0-03 
'^•30-0. 03 


4.8t-:i 


2.8l?i 


2 9+2-0 

^■y-0.7 


6 


44+0-04 
'-'•**-0.04 


0-28t°:«| 


0.059l«:~i 


[100] 



PI 97+0.03 
^■^'-0.03 
209+0-017 

U.ZUy_g gig 

n 1SS+0019 

0.1»»_g g3g 

20+0-02 

'-l-^U_g g2 

253+0-014 

U.ZOO_g gig 

944+0.003 
U.Z44_g gg3 



1.2^ 
2.2^ 
2.4^ 



+0.2 
-0.2 

r, + 0.5 

-0.3 
^+0.5 
-0.3 
2 6+0-3 
, -9+0.23 
l-'^-0.17 
, co+O.ll 
1-OO-0.07 



1.61 

2.08 
2.28 
2.62 
2.322 
1.95 



+0.02 
0.02 

+0.05 
0.04 
+0.04 
0.04 
+0.03 
0.03 
+0.017 
0.022 
+0.02 
0.02 



f-O.OOS 
-0.007 
1+0.02 
-0.02 
1+0.02 
-0.03 

0.688lHi^ 

o.eioiHJJI 

3+0.008 
-0.007 



0.385: 
O.54I 
0.59"* 



0.489: 



0.45 
0.87 
0.85 
0.83 
0.75 
0.89 



462.8/458 
242.9/217 
274.1/260 
559.2/563 
803.6/814 
690.7/656 



Note. — Col. (1): Observation index corresponding to Column (1) of Table [Tl Col. (2): Absorption column density in units of 102l cm 2. 
Col. (3): r is the photon index of the powerlaw model; r is the optical depth of the comptt model. Col. (4): A^pL is the normalization of 
the powerlaw model at 1 keV in units of 10~4 photons cm~2 s~l; Ac is the normalization of the comptt model in units of 10~4. Col. (5); 
Plasma temperature in units of keV; values in box brackets are fixed in the fitting. Col. (6): Inner disk temperature Tj^ of the diskbb model or 
seed photon temperature To of the comptt model in units of keV. For comptt + diskbb, To is set equal to Tin. Col. (7): i?in is the inner disk 
radius in units of 103 km; i is the disk inclination angle. Col. (8): Absorbed flux in 0.3-10 keV in units of 10~12 ergs~i cm~2 s~i. Col. (9): 
Unabsorbed luminosity in 0.3-10 keV in units of 1040 ergs s~i. Col. (10): bolometric luminosity (integration of physical models in 0.01—100 
keV) in units of 1040 ergs s~l. Col. (11): Best-fit E'nd degrees of freedom. All errors arc at 1 cr level. 



ton temperature Tq are equalized in the fits, i.e., assum- 
ing that the disk emission provides source photons for 
Comptonization. In some observations, the spectra show 
a featureless power-law tail at energies above 2 or 3 keV, 
leading to difficulty in finding the electron temperature 
Te, which is determined by the position of the high en- 
ergy bend. In those cases, we fixed at 100 keV to 
create a featureless power-law tail below 10 keV; higher 
or slightly lower temperatures do not affect other spectral 
parameters and the goodness of fit. 

Diagrams between pairs of parameters of interest are 
plotted in Figure [H The model being used to derive the 
parameters of that plot is indicated on the upper left 
corner of the panel (not including the absorption com- 
ponent). A tight, linear correlation was found between 
log Lx and the power-law photon index F when fitting 
the data with an absorbed power-law model. The best- 



fit relation is logLx oc (0.67 ± 0.03)F. The correlation 
coefficient between logLx and F is 0.995 with a chance 
probability of3.3xl0~^. The absorption column density 
Nn is also correlated with F, and seems to saturate when 
F > 2.6. 

When an additional disk component is added to ac- 
count for soft excesses, a correlation appears between 
log Lx and log Tin, with a coefficient of 0.95 and a chance 
probability of 3.5 x 10~^. Fitting with a power-law re- 
lation, we derive Lx oc jj^.iio.s ^pjiig resembles the in- 
trinsic evolution expected for a disk blackbody, except 
that the exponent of the relation is lower than 4. Tjn 
derived from the model diskbb is the color temperature, 
while the L oc relation is actually expected between 
the bolometric luminosit y and effective temperature Toff 
(jMcClintock et all l2008i ). We thus calculated harden- 
ing correction / = Tij^/Tcs as a function of luminosity 
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Fig. 1. — Diagrams between best-fit spectral parameters of NGC 5204 X-1 derived from different models (see Table[2)|- TLff is calculated 
from Tin with hardening correction (see texts for details). Dashed lines indicate best-fit relations; log(Lx) (0.67it 0.03)r, Lx oc T^^^^'''^ , 
and Lx '^cS^^^- "^^^ correlation coefficient r and chance probability p are: r = 0.995, p = 3.3 X IQ-^ between logLx and T, r 
p = 3.5 X fO"'' between logLx and logTin, and r = 0.92, p = 9.8 X fO"'' between logLx and logT^ff . 
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Continued. 



in Eddington units by compari ng the disk a tmosphere 
model bhspec (Davis et a.1. 200^ Davis fc Hu benv 2006) 
with the kerrbb model (jLi et alj 20051 ) assuming a non- 
spinning 100 Mq black hole (fitting with kerrbb to spec- 
tra generated based on bhspec), and then computed Toff 
assuming the highest observed luminosity is the Edding- 
ton limit. The computed / goes almost linearly from 
1.5 to 2.1 in the luminosity r ange o f 0.1 to 1 Eddington 
hmit (see iMcClintock et al.l ()2008l ) for a similar plot). 
After applying the hardening correction, we obtained 
Lx oc T^ff . The exponent of the relation between Lx 
and Tcff decreases if the hardening correction is computed 
at a lower Eddington level, e.g., we find an exponent of 
2.7 ± 0.9 if assuming the highest observed luminosity is 
half of the Eddington limit. This hardening correction is 
model dependent, but the basic behavior of the luminos- 
ity versus temperature is similar to that found for stellar- 
mass black hole X-ray binaries. The disk inner radius 
derived from the diskbb model is about 1000 km, with 
a tendency to decrease at high temperatures. A corre- 
lation also exists between fx and Tin, which strengthens 
confidence in the relation between Lx and Tin. A^h goes 
higher with Tin and Lx, indicative of a denser environ- 
ment when the accretion rate increases. The power-law 
component seems to be related to the thermal compo- 
nent, shown as a correlation between F and Tin. The disk 
emission contributes 0.11, 0.22, 0.19, 0.15, 0.23, and 0.23 
of the total luminosity in the 0.3-10 keV band, respec- 



tively for the six observations with a disk plus power-law 
model. 

The cool disk component peaks at energies where the 
absorption takes significant effect. Therefore, how to 
model the absorption must be considered carefully. We 
tried tbabs instead of wabs in XSPEC but no signifi- 
cant changes were found. We then tested ef fects caused 
by different abundances. IWinter et a.l ()2007f ) found that 
most ULXs including this one have near-solar abun- 
dances by measuring ox y gen a bsorption edge in XMM 
spectra. [Roberts et al.l ()2006l) found sub-solar abun- 
dance could marginally improve the fits to Chandra data 
of this s ource. The source is also found near young star 
clusters ()Goad et al.ll2002D . These results suggest that 
the abundance is close to or, perhaps, somewhat lower 
than solar. To see the effect of abundance, we modeled 
absorption within the Milky Way using wabs with A'h 
fixed at the Galactic value and the abundance set to solar 
and absorption within the host galaxy using vphabs with 
free A^h and the abundance set to 0.1 solar. The model 
consists of an accretion disk and a power-law. With sub- 
solar abundance, the derived luminosity is lower and the 
disk temperature is higher than with solar abundance 
(see Figure [5]). Fitting with a power-law relation, we ob- 
tain Lx oc Tjjj'^''^*''^, and a correlation coefficient of 0.97 
between log Lx and log Tin with a chance probability of 
1.1 X 10~^. Other abundances from 0.1 to 1 solar value 
will produce a slope in between. Therefore, the correla- 
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Fig. 2.— Lx versus T'm inferred from the diskbb plus powerlaw 
model with absorption in 1 and 0.1 solar abundances in the host 
galaxy, respectively. Data points for 1 solar abundance are the 
same as in Figure[T] With 0.1 solar abundance, the best-fit relation 
is L oc Tj^'^^" '^, and the correlation coefficient between log Lx and 
logTin is 0.97 with a chance probability of 1.1 X 10~^. 

tion between log Lx and log Tin appears robust even al- 
lowing for the uncertainty in the abundance of the host 
galaxy. 

If we model the spectra using an absorbed Comp- 
tonization model, the derived absorption column den- 
sity is rather low, often pegged at the Galactic value, 
which means no extra absorption is found in the host 
galaxy or internal to the binary system. This is some- 
what unphysical and hard to explain. In observations 1 
and 4, the best-fit Te is lower than 3 keV, as a result 
of spectral curvature. The optical depth in these two 
observations is also high, indicative of a cool, optically 
thick electron corona dominating the emission. However, 
the source in these two observations did not show dif- 
ferences on flux and absorption from the other four, in 
which it had power-law spectral tails implying emission 
from a hot, optically thin corona. The luminosity de- 
rived from the absorbed Comptonization model presents 
a positive correlation with the seed photon temperature, 
with a relation of Lx oc tJ-27±o.09 ^^^^ ^ yj.43±o.i5 

Then, we added a disk component in order to account for 
non Comptonized disk photons. With this model, a cool 
corona at low significance is also shown in observation 
5; if we fix Tg at 100 keV, the only increases 1.7 and 
no significant changes are observed on TVh and Tq, but 
the optical depth r becomes 0.11 due to the degeneracy 
between r and Tg. The luminosity, Lx, Lboi or Ldisk, 
becomes small when Tin goes up, which is obviously in 
conflict with the L oc T** relation. 

The disk blackbody plus power-law model is preferred 
for NGC 5204 X-1, because a correlation with moderate 
confidence level was found and a self-consistent interpre- 
tation can be made. However, replacing the power-law 
component with a physical (Comptonization) model fails 
to produce the anticipated correlation. This is probably 
because that the corona hides the innermost part of the 
disk when it is optically thick, and Tin = Tq is no longer 
valid in this case. Disconnecting Tin and Tq would lead 
to degeneracy in the model fits, making it difficult or 
impossible to find a unique global minimum, especially 
since the Comptonization model alone provides adequate 
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Fig. 3.— XMMPN spectra of NGC 5204 X-1 in observations 4 
(red) and 1 (blue). It is clearly seen that the high state spectrum 
has excesses over the low one at energies up to 3 keV, where ab- 
sorption takes very little effect. The low state spectrum cannot be 
modeled as the spectrum from the high state with extra absorp- 
tion. In the energy band above 3 keV, the two spectra have almost 
the same flux level, indicating the flux variation is due to a soft 
emission component. 



TABLE 3 

XMM OBSERVATIONS OF HOLMBERG II X-1 



Good Exposures (ks) 



No. 


ObsID 


Date 


PN 


MOSl 


M0S2 


1 


0112520601 


2002-04-10 


4.6 


10.1 


10.1 


2 


0112520701 


2002-04-16 


3.5 


5.3 


5.5 


3 


0112520901 


2002-09-18 


4.2 


6.5 


6.5 


4 


0200470101 


2004-04-15 


35.0 






Note. 


— Good Exposures are 


effective 


exposures after 



background screening. MOS data are not used for observa- 
tion 4. 

fits. A real physical model that includes the disk emission 
and Comptonization self consistently could be evaluated 
when hard X-ray data are available and tight constraints 
on the Comptonization parameters can be obtained. 

In Figure [3l two background subtracted energy spectra 
were plotted together, one from observation 4 with the 
highest flux among all, and the other from observation 1 
with the lowest flux. The two spectra were normalized in 
units of counts cm~^ s~^ keV~^ to account for different 
exposures and extraction areas. It is clear that the high 
flux spectrum is above the low one at energies below 
3 keV, and they are almost at same levels above 3 keV. 
This means that the flux variation is from a soft emission 
component, and is not caused by absorption that takes 
almost no effect at energy around 2-3 keV with A'h ^ 
10^1 cm-2. 

2.2. Holmberg II X-1 

Holmb erg II X-1 is surrouncied by optical and radio 
nebulae (iP akull k MiriMfi|[200l fMilier et al.|[2005[ ). The 

optical nebula has a roughly isotropic morphology, and is 
thought to be photoionized r equiring an input X -ray flux 
of at least 4 x 10^^ ergs s'^ (|Kaaret et al.ll2004D . These 
indicate that it is a truly ULX with little beaming. The 
source has been observe d four times wit h XMM (Table [31 
iDe wangan et a l. 2004; G oad et "ani2006[ ). and three times 
with Chandra including one 100 ks observation (Ob- 
sID 5933) which has never been reported. Here, we re- 
analyzed all archival XMM data, and included the 100 ks 
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TABLE 4 

Best-fit parameters of Holmberg II X-1 























No. 

(1) 


(2) 


r/r 

(3) 


iVpL/Afc 
(4) 


To Tin/To 

(5) (6) 


Rin V COS i 

(7) 


fx 
(8) 


Lx 
(9) 


(10) 


(11) 


Model: wabs*powerlaw 


1 
2 
3 

4 


-0.04 
-1 OO+0.04 
'-■J-i-0,04, 
1 45+0.08 

-, Ro+0.02 


9 70+0-02 
^- '"-0.02 
9 47+0.03 
^■^'-0.03 
05+0.06 
■J-'-' -^-0.05 
9 700+0.012 


27 4+0-4 
21 0+0-4 
D.ci_Q 2 

27.9l0;2 






6 56+0-05 
"■''"-0.05 

6 20+0-0'' 
'^•^'-'-0.07 
1 49+0.03 
l-+^-0.03 
6 6S+0-03 


1 00+0. 03 
l-°'-'-0.03 

1 40+0-02 
l-+'-'-0.02 

46+0-02 
^•^"'-0.02 
1 Sfifi+0-018 




682.2/653 
560.3/508 
220.7/196 
817.7/700 










Model: wabs(powerlaw 


+ diskbb) 










1 
2 
3 
4 


1 41+0.07 

-1 90+0.09 

1 3Q+0.14 
^■'-'^-0.12 
1 54+0.04 


9 47+0.05 
^•^'-0.09 
9 90+0.05 
^•^■^-0.05 
9 oA+0.09 

9 55+0.02 
^•^'^-0.03 


20.2+2 4 
16 

^°-^-i.o 

5 9+0.5 
23 1+0-'' 


U.ZaZ_Q gjg 

n 995+0.019 
u.zzo_(, Q20 

17+0-02 

^•l'-0.02 
Q 260+0-018 


1 29+0-10 

2 4+0-'' 
^•^-0.4 

2 5+1-3 


6 67+0-06 
D-D'-0.06 

"•^'-'-0.08 
1 45+0-03 

6 74+0-03 
'^-0.03 


55+0.07 
1-33-0.07 
-I 00+0. 05 
^•'^*-0.05 
43+0.04 
"•^■^-0.04 
-I 71+0.04 
1- '1-0.04 




656.8/651 
518.0/506 
211.5/194 
756.4/698 


Model: wabs*comptt 


1 
2 

3 
4 


U.D'_o.09 

n 76+0-08 

^•'"-0.12 
'J- '^-0.18 

n •50+0.06 


077+0.003 
'-0.003 
,07+0.008 

OOP.+0.004 

U.UOO_g QQ3 

4 •5+0-3 


fifi+0-10 

'^•"0-0. 06 
fiq+0-08 
"•"^-0.09 

90+0.09 

20.11 


[100] 0.139l°;««« 

[100] o.iwtlZl 

[100] 0.103lO;°i;5 
2q+0-4 n 155+0.005 
^•='-0.3 "■^■^"'-0.005 




6 fio+0.05 

90+0.10 

0- ^*-0.05 

1 44+0-03 

1- ^^-0.02 
70+0.03 
"• ''-'-0.03 


-1 19+0.06 
^•1^-0. 04 
1 09+0-03 

i-uy-0.05 

29+0-03 
^•^^-0.03 

1 0,5+0.020 

±.UCiO_o 019 


1.35 
1.41 
0.34 
1.12 


643.3/652 
524.1/507 
207.8/195 
702.7/698 










Model: wabs(comptt + 


diskbb) 










1 
2 

3 
4 


n SI +0.10 

U.8i_o.l4 

n qfi+0 09 

'^•^'-0.13 

n 75+0.05 


n 076+0-003 

U.UrD_Q QQ3 

5 7+1.0 
' -2.0 
7 q+3-2 
'•='-1.3 
4 fj+0.5 


0.54t»;J^ 

13.1^; 

9+1 


[100] 0.159«;«^; 
2 4+11 -,0^+0.020 

1 1+0-4 ,ot-+0.027 
1-1-0.2 '-'-1°''-0.018 

2 6+0-4 0,09+0.011 


4 1+0-' 
+-1-0.9 

4 2+1-0 
+■^-0.7 
q+1-3 

3 6+0-2 


6 66+0-05 

6 ■51+0-09 

0- '^l-O.lO 
1 41+0.03 

1- *l-0.03 
79+0.03 
"•'^-0.03 


1 17+0.05 
l-l'-0.05 
1 17+0.05 
l-l'-0.04 

O1+0.03 
^^•■^1-0.02 

1 10+0.02 
l-l°-0.03 


1.53 
1.51 
0.46 
1.43 


641.2/651 
514.5/505 
203.2/193 
677.6/697 



Note. — Columns are the same as Table [2] 

Chandra data for timing analy sis. We adopted a distance 
of 3.39 Mpc to the host galaxy (jKarachentsev et al.ll2002f ) 
and an absorption column density of 0.365 x 10^^ cm~ 
in our Galaxy. 

The radius of the circular region used to extract source 
spectra is 29" for the first two observations and 30" 
for the last two, respectively . The 4th observation has 
been reported in i Goad et al.l (|2006'). in which the spec- 
tral analysis was made using a combination of 0.3-6 keV 
MOS data and 0.7-10 keV PN data in order to avoid cal- 
ibration issues at other energy ranges. Here, we simply 
used 0.3-10 keV PN data for spectral fitting, since the 
st atistics are go od en ough. 

iGoad et all ()2006[ ) suggested that the absorbing 
medium in the host galaxy has an oxygen abundance 
of 0.56 solar from spectral fits to RGS data. However, 
IWinter et a.l (pool found near-solar abundance of oxy- 
gen from EPIC data. As we have discussed previously, 
lower abundances systematically change the inferred lu- 
minosity and disk temperature but will not affect the 
confidence level of the correlation. We therefore use so- 
lar abundance in the fits. 

The same four models used for NGC 5204 X-1 were ap- 
plied to Holmberg II X-1, and the best-fit parameters are 
listed in Table S) The power-law model provides poor fits 
at energies around 0.5 keV, which could not be improved 
by adding a disk component (see Figured]). In contrast, 
the Comptonization model, whether or not with a disk 
component, produces uniformly good fits in the full band. 

Diagrams between parameters of interest are plot- 
ted on Figure [5] Interestingly, the luminosity inferred 
from a disk plus power-law model goes higher when 
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Fig. 4. — Data to model residuals for the 4th observation of 
Holmberg II X-1 in units of cr. Different panels are for different 
models with names shown on the upper right corner. Models con- 
taining a power-law component are unable to fit the data around 
0.5 keV. Comptonization models provide good fits throughout the 
XMM energy band. 

the disk temperature becomes hotter, with a relation of 
Lx oc TJ^'^^^'^, and a correlation coefhcient of 0.84 with 
a chance probability of 0.16. This is consistent with the 
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Fig. 5. — Diagrams between best-fit spectral parameters of Holmberg II X-1 derived from different models (data points adopted from 
TableHJl. Dasfied lines indicate best-fit relations; Lx oc T^;^*^'^, and Ljisk oc 7^4 i±2-2 



intrinsic evolution of disk blackbody emission, i.e., lumi- 
nosity changes with the 4th power of the temperature. 
However, the four data points provide only weak evi- 
dence for the correlation. More observations are required 
to test it. 

The Comptonization model provides better fits than 
the power-law model, especially at energies around 
0.5 keV. For the first three observations, the electron 
temperature is not well determined due to the lack 
of curvature in the spectrum at high energies and is 
thus fixed at 100 keV. No significant correlation can 
be seen between luminosity and seed temperature for 
a single Comptonization model. When we include a 
disk component, observation 2 and 3 also present pos- 
sible high energy spectral curvature, and a plausible cor- 
relation between Ldisk and Tin- The best- fit relation 
is Ldisk oc T|^'^*^'^, and correlation coefficient between 
log Ldisk and log Tin is 0.93 with a chance probability of 
0.072. The best-fit inner radius of the disk is consistently 
around 4 x 10^ km. 

In Figure [SI we present lightcurves in the 1-4.5 keV 
energy range of Holmberg II X-1 from the longest Chan- 
dra (ObsID 5933) and XMM (no. 4) observation. For 
the Chandra data, the lightcurve is extracted from a 3cr 
elliptical region around the source found by wavdetect, 
and binned at a timescale of 800 times the frame time 
(about 2.03 ks). No background is subtracted for the 



Chandra data due to its small fraction. For the XMM 
data, the lightcurve was extracted from the PN data only, 
corrected for background using the same region for spec- 
tral analysis, and binned at a timescale of 10,000 times 
the frame time (about 0.73 ks) at intervals without tim- 
ing gaps. Strong variations are obviously seen from the 
lightcurves. Periodograms were calculated to examine 
typical timescales of the variation. A relatively narrow 
peak around 20 ks is pronounced, corresponding to mul- 
tiple spikes in lightcurves at the same timescale. Due 
to the limited duration, we cannot determine whether or 
not the oscillation is a periodic signal or possibly a red 
noise fluctuation. We also examined variability at other 
energy ranges, and found they are not as prominent as 
in the 1-4.5 keV band. 

2.3. IC 342 X-1 and X-2 

IC 342 X-1 and X-2 are famous because they are 
the first two ULXs in which state transitions were dis- 
covered (jKubota et al.l l2001h . and are thus ideal ob- 
jects to st udy here. IC 342 X-1 resides in an opti- 
cal nebula (iPakuUfc Mirionill2002t [Roberts et all 120031 : 
iFeng fc Kaaretl[2008h . which is similar to Holmberg H X- 
1. The galaxy has been observed four times with XMM 
(Table[5|) , and the last three observations have never been 
reported. 

The radius used to extract spectra for X-1 is 25" for the 
1st observation and 30" for the other three; for X-2, it is 
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Fig. 6. — 1-4.5 keV lightcurves and periodograms of Holmberg II X-1 obtained from Chandra ACIS and XMM PN observations. Timing 
gaps have been removed and background has been extracted for XMM PN data. 



TABLE 5 

XMM OBSERVATIONS OF IC 342 









Good Exposures (ks) 


No. 


ObsID 


Date 


PN 


MOSl 


M0S2 


1 


0093640901 


2001-02-11 


4.8 


9.5 


9.5 


2 


0206890101 


2004-02-20 


10.3 


20.4 


21.4 


3 


0206890201 


2004-08-17 


17.1 


23.4 


23.4 


4 


0206890401 


2005-02-10 


4.6 


7.3 


7.6 



Note. — Good Exposures are effective exposures after 
background screening. 



20" for the 2nd observation and 32" for the rest. For X-2 
in observation 1, the MOS data are not used because the 
source is located right across the CC D gap. We adopt a 
distance of 3.3 Mpc to the host galaxy (|Saha et al.ll2002D . 
and an absorption column density of 3.08 x 10^^ cni~^ in 
our Galaxy. 

We tried three different models to fit the spectra of 
X-1; the disk plus Comptonization model is not used be- 
cause a global minimum cannot be found for this model 
in some observations due to the weakness of the disk com- 
ponent. The best-fit parameters are listed in Table [51 
If we fit the data with an absorbed power-law model, 
the four observations appear in two groups: a low hard 
group and a high soft group (see Figure [7]). The absorp- 
tion column density is also correlated with the source 



luminosity. When the disk component is added, the disk 
temperature seems to be inversely scaled with the lumi- 
nosity with a best-fit relation of Lx k T;^^'^^°'^, which 
excludes a L oc T"* relation. Fitting with the Comp- 
tonization model, a cool, optically thick corona is derived 
as a consequence of spectral curvature at energies close 
to 10 keV. This is especially evident in the 2nd observa- 
tion, in which the Comptonization model significantly 
improved the fits relative to the disk plus power-law 
model. However, no correlation can be found between 
spectral parameters with the Comptonization model due 
to the relatively large errors. 

The shape of the energy spectrum of X-2 seems to be 
abnormal. We plotted the source spectra in the 2nd and 
3rd observations in Figure [H as well as residuals after 
fitting with an absorbed power-law model to the data in 
2-10 keV. Excesses at low energies are seen in the residu- 
al s. Unlike usual soft excess es found in ULXs (e.g.. Fig. 2 
in lGonqalves fc Sorial[2006h . these excesses peak around 
1 keV and have a minimum around 0.5 keV; they cannot 
be fitted by a cool disk component but can be explained 
by thermal plasma emission (apec in XSPEC) down to 
0.5 keV. Therefore, we first tried a power-law model and 
an accretion disk model modified by absorption to fit the 
data in the energy range of 1-10 keV to avoid the soft 
excesses. Since an apec model in addition to the power- 
law is able to fit the data down to 0.5 keV, we then 
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TABLE 6 
Best-fit parameters of IC 342 X-1 



No. 

(1) 


Nil 
(2) 


r/r 

(3) 


AfpL/ATc 

(4) 


To 
(5) 


Tin /To fiinVcOSi 

(6) (7) 


fx 
(8) 


Lx 
(9) 


-^'bol 

(10) 


xVdof 
(11) 


Model: wabs*powerlaw 


1 
2 
3 
4 


c- O+0.3 
3-«-0.3 

6.32l°;l^ 

8.41^2 


1 65+0-05 
^^-''0-0.05 

^•'J^-0.02 
1 oq+0.03 
^•''•-'-0.03 

1 89+°"^ 
^•'"'-0.03 


'^■'J-0.3 
ig 7+0.5 

"•='-0.2 
21 






9 CI +0.07 
^■^^-0.07 
c 0(5+0.05 
^•■3D_o.05 

2 64+°"'* 
^•°^-0.04 
7 no+0.10 
'•'^''-0.10 


Q 470+0.012 

u.'i(0_o 012 
1 41+0.02 
^•^^-0.02 

n 560+° °°^ 

U.ODU_o 009 

1 64+°-°'* 
-l-°*-0.03 




144.0/137 
750.2/667 
452.0/435 
421.0/374 


Model: wabs (powerlaw + diskbb) 


1 
2 
3 
4 


°-^-i.o 

ii.stj:? 


1 58+0-11 

i-a8_o.io 

9 17+0.05 
-0.03 
. 7g+0.05 
'"-0.05 

^■^^-0.04 


4-6li;^ 

24.6tl:l 
"-■^-0.5 

25.8j:?-0 




0.281°:^^ 0.9tg;^ 
„ -,30+0.007 94+13 

n 9R+0.05 1 9+1.0 
^•^°-0.04 ^-^-0.5 
n 197+0.010 ■31+24 
"•^^'-0.014 '-'^-10 


2 54+0-08 
^■O''_0.08 
t- 01 +0.05 
"•■^^-0.05 

2 67+° °'* 
^■°'-0.04 

6 g-s+o-io 


n 57+0-28 

4 1 + 1-8 
^■^-1.0 
56+0.10 

U-DD_0.06 

4 6+2" 




140.9/135 
734.4/665 
437.4/433 
409.8/372 


Model: wabs*comptt 


1 
2 
3 
4 


4 7+1-0 
^•'-0.8 
7 Q+0-2 

c 9+0.4 
'•^-0.3 


' -^-2.4 
7 Q+l.l 

5 Q+0.5 

7 9+0.6 
'•^-0.6 


o 7+1.4 
24+149 

4.01?:^ 


9 O+3.0 

9 14+0.17 
^■1^-0.13 
04+1.2 
■^•^-0.5 
9 9+0.3 
^•^-0.2 


n 90+0.05 

^^•^■^-0.07 
14+0-14 
^•1^-0.14 

r, 99+0.02 

^•^^-0.03 
U-UD_o 06 


9 47+0.11 
^■^'-0.09 
c: 91 +0.06 
"■^1-0.06 

2-64t°:°l 

6.80t°:- 


n 40+0-04 
1 n+0-12 

1-11_0.22 
n 455+0.019 
U.400_o 017 

^•^^-0.03 


0.51 
1.26 
0.58 
1.83 


140.6/135 
711.2/665 
439.8/433 
409.8/372 



Note. — Col. (1): Observation index corresponding to Column (1) of Table [5] Col. (2): Absorption column density in units of 
10^^ cm"'^. Col. (3): F is the photon index of the powerlaw model; r is the optical depth of the comptt model. Col. (4): A'^pl is the 
normalization of the powerlaw model at 1 keV in units of 10^^ photons cm"^ s~i; A^c is the normalization of the comptt model in units 
of 10~4. Col. (5): Plasma temperature in units of keV; values in box brackets are fixed in the fitting. Col. (6): Inner disk temperature Tj^ 
of the diskbb model or seed photon temperature To of the comptt model in units of keV. Col. (7): ifin is the inner disk radius in units of 
10^ km; i is the disk inclination angle. Col. (8); Absorbed flux in 0.3-10 keV in units of 10"*'^ ergs~^ cm"^ s"*. Col. (9): Unabsorbed 
luminosity in 0.3-10 keV in units of 10*° ergs s~^. Col. (10): bolometric luminosity (integration of physical models in 0.01-100 keV) in 
units of 10**° ergs s^^. Col. (11): Best-fit and degrees of freedom. All errors are at 1 cr level. 



TABLE 7 

Best-fit parameters of IC 342 X-2. The spectral fits were performed in the 1-10 keV range for the wabs*powerlaw and 
wabs*diskbb MODELS, BUT IN THE 0.5-10 keV range for THE wabs (powerlaw + apec) MODEL. 



No. 

(1) 


(2) 


r 

(3) 


iVpL 
(4) 


Tin/kT RiW COS i/N A 
(5) (6) 


fx 
(7) 


(8) 


i'bol 

(9) 


xVdof 
(10) 


Model: wabs*powerlaw 


1 
2 
3 
4 


9 0+0.3 
^•'-'-0.2 
9 co+0.07 
^•"•-'-0.06 

1 70+0.08 
^•'^-0.04 

2 14+°-2i 
^•1^-0.16 


1 79+0-14 
^- "'-0.13 
1 71+0.03 
^•'^-0.03 
1 07+0.04 
^•"'-0.04 

1 64+°-" 
^•°^-0.08 


.,+2.0 
8-3_i.6 

30.2j;i'^ 

Q+0.3 
0.3_Q 2 

6 4+*-2 
°-^-0.8 




r, c:o + 0.12 
^■0»_0.11 

10.55l°;i? 
2.50l°:°l 
2.59t°:°'8 


51+0-03 
"•"1-0.03 

9 ns+0-03 

^•U»-0.03 
422+°"°'' 
n 485+0-020 
0.485_o 018 




45.0/39 
717.2/605 
350.6/316 

74.5/84 


Model: wabs*diskbb 


1 
2 
3 
4 


1 5+0.2 

1 09+0.04 
l-»^_0.04 
1 ocr+0.08 
-■-■■J ''-0.04 

1 42+°-l2 
-'-•*'^-0.15 






917+0.22 n+O-6 

■^■'-'-0.18 ■^■'-'-0.5 
9 50+0-06 4 54+0.19 
^-OU_o 05 "i-o^-o.ig 

04+0.11 1 01+0.13 

■^"^^-0.19 J^-^1-0.07 
r,p,+0.27 91+0.2 
^•°1-0.15 ^-1-0.3 


2.4lt°:i? 

io.i4i°:ii 

9 40+0.04 
^■^"-0.04 
9 40+0.09 
^■^^-0.09 


0.408t°;°}^ 

l.707t°;°l? 
n 070+0. 006 

u.olo_o 006 
OQ7+0.012 

u.jy^_o 012 


0.50 
2.19 
0.55 
0.52 


34.9/37 
591.7/596 
303.6/311 

70.8/83 


Model: wabs (powerlaw + apec) 


2 
3 


99+0.10 
•""■^^-O.IO 
9 c:o+0.15 
^-3»-0.15 


1 01 +0.04 
l-yi-0.04 
1 ™+0.05 
l-"^-0.06 


44+1 
5 Q+0-6 


009+0-008 ^1+33 
U.Uyz_o 003 "1-27 

126+°"l'^ 4+°" 
U.iZD_ooi2 'J-^-0.3 


10.4lt°- 
2.47l°:°l 


3 6+°''' 
"■O-0.3 

60+°"'' 
"-""-0.05 




603.2/608 
327.6/320 



Note. — Col. (1): Observation index corresponding to Column (1) of Table[5] Col. (2): Absorption column density in units of 10^^ cm 
Col. (3): Photon index of the powerlaw model. Col. (4): Normalization of the powerlaw model at 1 keV in units of 10~4 photons cm"^ s~i. 
Col. (5): Tjji is the inner disk temperature of the diskbb model; kT is the plasma temperature of the apec model. Both are in units of keV. 
Col. (6): iJin is the inner disk radius in units of 10^ km and i is the disk inclination angle for the diskbb model; A'^a is the normalization 
of the apec model. Col. (7): Absorbed flux in 1—10 keV in units of 10~^^ ergs~i cm"^ s~i. Col. (8): Unabsorbed luminosity in 1—10 
keV in units of 10*° ergs s~i. Col. (9): bolometric luminosity (integration of physical models in 0.01-100 keV) in units of 10*° ergs s~i. 
Col. (10): Best-fit and degrees of freedom. All errors are at 1 cr level. 
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Fig. 7. — Lx versus Tin of IC 342 X-1 derived from the disk 
plus power-law model (data points adopted from Table IBJ . The 
four data points represent a relation of Lx °^ jj-2-8±0.7 




2.0 



2.5 3.0 



3.5 



4.0 



Fig. 8. — Lboi versus Tin of IC 342 X-2 derived from the absorbed 
disk model (data points adopted from Table (Tj. The spectral fits 
were performed in the energy range of 1-10 keV to avoid the soft 
excess whose origin is unclear. 



applied such a model to the 2nd and 3rd observations 
in the 0.5-10 keV range; the 1st and 4th observations 
do not have enough counts to allow significant detection 
of the apec component. We note that we cannot find 
a physical model that is able to adequately fit the data 
down to 0.3 keV for the 2nd observation. Below 0.5 keV, 




0.5 1 2 

Energy (keV) 

Fig. 9.— XMM PN spectra of IC 342 X-2 in observations 2 
(blue) and 3 (red). Measured spectra are shown in the top panel, 
and their residuals (in units of a) after fitting with a power-law 
model in 2-10 keV are shown in the middle and bottom panels. An 
soft excess between 0.5-1.5 keV is obvious in both spectra, which 
can be adequately fitted by the thermal plasma model apec. The 
high state spectrum, however, present flux excess below 0.5 keV, 
which seems to be a super soft emission component and not subject 
to the interstellar absorption. 

there seems to be another emission component that is 
not subject to the absorption. All best-fit spectral pa- 
rameters for different models are listed in Table [71 To be 
consistent, fx and Lx are all quoted in the energy range 
of 1-10 keV. 

Two distinct flux states are found in the four observa- 
tions, but no correlation is shown between spec t ral pa - 
rameters. Unlike that reported in iKubota et al.l (|2001h . 
the low state and the high state has almost the same 
spectral shape except at energies below 0.7 kcV (Fig- 
ure [21) . Fitting with a single accretion disk model indi- 
cates that the source stays in the thermal dominant state, 
in which the L oc T'^ relation should robustly exist. How- 
ever, the derived Xboi and Tin strongly violate the L (x 
relation (see Figure ^ ; fitting with a L oc function 
results in = 120.9 with 3 degrees of freedom. This 
is solid evidence against the hot disk model, indicating 
the source emission is not from a standard, hot accre- 
tion disk onto a stellar-mass black hole. The apec model 
is able to explain the soft excess between 0.5-1.4 keV. 
In the 1-10 keV range, the apec component contributes 
38% and 22% of the total luminosity, respectively for the 
2nd and 3rd observation. 

2.4. Antennae X-11, X-16, X-42, and X-44 

The Antennae (NGC 4038/4039) are a pair of merging 
galaxies, which have a high star formation rate and har- 
bor a wealth of br ight X-ray sources as well as ULXs 
(jZezas et al.l l2002a( ) . With archival Chandra observa- 
tions, four ULXs, X-11, X-16, X-42, and X-44, are 
found to have particula rly high luminosity and variability 
(jFeng fc Kaaretll2006bO . We proposed six XMM obser- 
vations of the Antennae (PI: H. Feng) in order to in- 
vestigate spectral variabilities of these UL Xs. Another 
two ar chival observations are also included ([Miller et al.l 
l2004bH . 



12 



Feng and Kaaret 



TABLE 8 

XMM OBSERVATIONS OF THE ANTENNAE 









Good Exposures (ks) 


No. 


ObsID 


Date 


PiN 


MOSl 


MUS2 


1 


0085220101 


2002-01-18 


14.9 


19.9 


20.0 


2 


0085220201 


2002-01-08 


33.6 


51.2 


51.4 


3 


0500070201 


2007-06-09 


9.1 


13.5 


13.6 


4 


0500070301 


2007-06-24 


19.2 


24.1 


23.9 


5 


0500070401 


2007-12-10 


8.4 


18.4 


19.7 


6 


0500070501 


2007-07-09 


14.6 


24.2 


10.3 


7 


0500070601 


2007-12-18 


17.4 


25.1 


24.4 


8 


0500070701 


2007-12-26 


21.7 


31.1 


31.0 



Note. — Good Exposures are effective exposures after 
background screening. 



Due to the crowded field and strong background emis- 
sion from the interstellar medium (ISM), we extracted 
spectra from smaller regions than usual, which are 12" 
radius for X-11 and X-16 and 10" radius for X-42 and X- 
44. As mentioned above, the small extraction region does 
not affect the inferred flux. The PN data in observations 
1, 2, 5, 7, and 8 are not used for X-11, because it hes on 
the CCD gap. No background subtractions because the 
small extraction region results in negligible instrumen- 
tal background. The ISM does contribute a background, 
but this varies with position in the galaxy. To account 
for the ISM emission, we model its spectrum using the 
thermal plasma model apec with parameters adjusted for 
each source. The distance to the Antennae is controver- 
sial: ~20 Mpc from recession velocity, 13.3 Mpc derived 
by measuring the t ip of the red giant branch (TRGB; 
ISaviane et al.ll200"8[ ). and 22 Mpc based on a type la su- 
perno va and re-analysis of the TGRB (iSchweizer et al.l 
l2008f ). In our previous paper (iFeng fc Kaaretll2006bn . 
we adopted a distance of 13.8 M pc, which is estimated 
from TRGB (jSaviane et al.ll2004l ). Here, we adopt the 
most recent result of 22 Mpc, which leads to higher lu- 
minosities by a factor of 2.5 than before. The absorption 
column density in our Galaxy is 0.314 x 10^^ cm~^. 

The ISM emission for each source was determined it- 
eratively. First, we fit every spectrum using a power-law 
plus apec model with all parameters free except the apec 
abundance (set to solar) and redshift (set to zero). Then, 
we fixed the plasma temperature for each source at its 
mean value averaged from different observations and fit 
the spectra again to find new normalizations. Finally, we 
fixed the apec plasma temperature to the mean temper- 
ature and the normalization to the mean normalization. 
The resulting apec model for each source served as a 
background spectrum for that source. We note that the 
final mean values are consistent with the best-fit values at 
the first step for each single observation, which is reason- 
able as a result of steady ISM emission. The temperature 
is 0.43, 0.38, 0.44, and 0.45 keV, and the normalization is 
9.6, 5.6, 13, and 3.2 (in units of 10~^ times the apec nor- 
malization; sec XSPEC manual), respectively for X-11, 
X-16, X-42, and X-44. According to the quahty of the 
data and small contribution of the ISM emission, slight 
changes of those values will not affect the best-fit values 
of other parameters. 

We tried a single power-law model subject to interstel- 
lar absorption to fit the data. All best-fit parameters 
are listed in Table [9] We plotted Lx versus power-law 



photon index F in Figure [10] for the four sources. It 
is interestingly shown that all sources, especially X-11, 
X-16, and X-42, display a constant F despite dramatic 
changes of the luminosity. X-44 shows some trends of 
hardening when the source becomes bright, which is con- 
sistent with results from previous Chandra o bservations 
(|Fabbiano et all 120031: iFeng fc Kaaretll2006bf) . X-16 ex- 
hibits an extremely hard spectrum with F = 1.1 ~ 1.4, 
and a low absorption column density close to the Galactic 
value. 

3. DISCUSSION 

In the past decade with Rossi X-ray Timing Explorer 
observations, an observational description of spectral 
evolution in black hole binaries has been established 
with the identification of four states (quiescent, hard, 
thermal dominant, and steep p ower-law) and transitions 
between them (for reviews see iMcClintock fc RemillardI 
l200^lR,emillard fc McClintockll200(f . Each state is char- 
acterized by its unique spectral, timing, and multiwave- 
length properties, and linked to a specific physical con- 
figuration of the accretion flow. A standard, phenomeno- 
logical spectral model consisting of a disk blackbody plus 
power-law is widely adopted to describe the X-ray spec- 
trum of black hole binaries. Here, we discuss our results 
on ULXs by comparing them with Galactic black hole 
binaries in terms of these spectral states. 

3.1. The Cool Disk Model 

The correlation between the luminosity and disk tem- 
perature of NGC 5204 X-1 suggests that the soft ex- 
cess in the spectrum can be interpreted as emission 
from an accretion disk. The L cx relation found 
between the luminosity and inner disk temperature is 
similar to that seen at the highest luminosities ob- 
served from_JCTE_J1550^64 in its 1998-1999 out- 
burst (jKubota fc Dond l2004^ and from H1 743-322 in 
its 2003 outburst (jMcClintock et al.l l2007f ). The ef- 
fective temperature observed from an accretion disk 
is modified by 'spectral hardening' due to the disk 
atmosp here. Using stat e -of-th e-art disk atmosphere 
models, IMcClintock et all ()2008D analyzed the data for 
XTE J1550-564 and H1743-322 and showed that the 
deviations from the L oc T'' seen at high luminosity when 
the luminosity is plotted against the observed disk tem- 
perature are largely removed when the spectral harden- 
ing is taken into account and the luminosity is plotted 
against the effective disk temperature. Our analysis pre- 
sented above, shows that application of a spectral hard- 
ening correction to the data for NGC 5204 X-1 produces 
a similar effect: the luminosity versus temperature re- 
lation for the effective disk temperature is close to the 
L oc relation expected for an accretion disk. There- 
fore, we interpret the source emission of NGC 5204 X-1 as 
due to a standard or quasi-standard accretion disk with 
the observed spectrum modified by spectral hardening. 

The mass of the compact object can be derived un- 
der this scenario. Assuming the source peak luminosity, 
^ lO^*' ergs s~^, is under the Eddington limit, which is 
1.3 X 10'^* ergs s~^ for a solar mass black hole, a lower 
limit on the mass can be calculated as M > 80 Mq. The 
inner radius derived from the disk model has a maxi- 
mum of about 1000 km and slowly declines at high tem- 
peratures to about 600 km. Similar to XTE J1550— 564 
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TABLE 9 

Best-fit parameters of Antennae X-11/16/42/44 with an absorbed power-law model, wabs(powerlaw + apec), where apec 

PARAMETERS ARE FIXED TO SERVE AS A BACKGROUND SPECTRUM. 



No. 


Nn 


r 


JVpL 


/x 


ix 


xVdof 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


X-11 



1 


0.31<0-^i 


1 61+°-" 

'-■"'--0.09 


n sn5+o-042 


9,0+0.014 


1 oo+O.lO 
-^■'-'•^-0.09 


20.4/17 


2 


l-J^D-O.ll 


1 76+0-05 
1- '"-0.05 


'-'■^''-0.04 


301+° °" 


9 11+0-07 
^■^^-0.07 


67.1/67 


3 


1 24+0-1'' 


1 'ifi+o-os 

^■^"-0.08 


39+0-04 


262+°-°" 


1 76+0-09 

1- 'D-0.09 


23.7/29 


4 


n 75+0-10 

"•'°-0.11 


1 fiS+0-06 


34+0-02 


2oq+°°i° 


1 Qo+0.05 
^■3»_0.05 


56.2/59 


5 


0.40<l-95 


1 n+0.3 

^•■^-0.3 


Q ,3+0.08 
"■^'-'-0.03 


13+0-02 
'^■^•^-0.03 


79+0-11 
'^-0.10 


2.4/6 


6 


1 24+0-i« 


1 79+0.13 
^- '^-0.13 


n 90+0.04 

"■^^-0.03 


147+0.013 
'^•^^'-0.012 


1 (,^+0.06 


39.6/29 


7 


32<0-66 


^■^'-'-0.13 


n 150+0-029 


134+0.013 

U.X01_g gj^4 


r, 01 +0.07 
U.»i_g.g7 


19.1/12 


8 


0.33<0-5'3 


-I oQ+0.14 
l.»3_o 


n 1H-3+0.036 
U.i»3_Q 


n 105+0-010 

U.lUO_g ggg 


66+° °^ 
"■""-0.05 


29.8/16 



X-16 



1 


n 47+012 
'J-^'-o.ii 


1 26+°-°° 
l-^"-0.07 


n 97+0.02 

^^•^'-0.02 


0.284i°;°l^ 


1 79 + 0.08 
1- '^-0.08 


31.9/44 


2 


n 46+°-°'' 

"■^"-0.07 


1 oc:+0.03 
1-JO_0.04 


0.318l°:°;| 


0.298t°°°^ 


-, 09+0.05 
^■°^-0.05 


120.4/119 


3 


0.3l<0.38 


1 90+0.14 

i-^'-'-o.io 


0.133t°:°i^ 


n -1^:9+0.015 

U.iOy_g gig 




13.8/12 


4 


0.3l<0.39 


1 99+0.09 
1-^^-0.09 


0.0981°:™? 


n 11 9+0.008 
U.J-J-^_g gig 


n fis+0-04 

U.b»_o g4 


14.4/18 


5 


0.3i<0.50 


1 11+0.08 
1-11-0.05 


n 910+0.022 
u.^ie_o gii 


0.283l°°l? 


1 fiS+0-09 
1-D»_0.09 


24.4/29 


6 


0.86l°:it 


1 01 +0.07 
I'-'l-O.OS 


r, 99 + 0.02 

^■^^-0.02 


0.215l°°ll 


1 34+0.07 
l"5^-0.07 


16.6/30 


7 


0.78l°:il 


-, 99+0.05 
1-^^-0.05 


30+0-02 
"■•^"-0.02 


0.30ll°°ll 


-, 07+0.07 
l-»'_0.07 


67.2/57 


8 


Q 37<0.44 


1 90+0.05 
l-^'-'-0.05 


242+°-°l'5 

U.Z4Z_g gj3 


0.280t°;°ll 


1 co+0.06 

i-o^-o.oe 


73.0/68 



X-42 



1 




1 7,,+0.10 

1- '•^-o.ii 


25+°-°^ 
'^•^'^-0.03 


144+0-012 


96+0-05 
"■^"-0.05 


31.0/31 


2 


0.63l°:°l 


1 66+°°'' 

l-DD_g gg 


241+0-017 

U.^41_g gj^g 


n -|t;5+0.006 

U.100_g ggg 


1 nn+0.03 
l-'J^-0.03 


83.4/86 


3 


1 4+0-3 

1-^-0.7 


9 9+0.3 
^•^-0.5 


90+0. 06 
"■^'^-0.09 


0.083t°°lt 


65+°-l'' 
1-bO-0.07 


14.7/12 


4 




9 91 +0.09 
^-^1-0.11 


r, 90+0.03 

"-^'-'-0.03 


0.08ll°°°i 


67+° °'' 
1-b' -0.06 


46.6/31 


5 


o.65i°:iI 


1 70+0-10 
1- "J-o.io 


27+° °'' 
"■^' -0.03 


0.168l°:°l? 


1 10+O.O6 
1-1^-0.06 


35.3/26 


6 


0.53l°:i? 


-I 7C+0.18 

1- '"'-0.17 


19+0-03 
U.ia_g.g3 


„ 117+0.011 
"•ii'-o.oio 


n 7cr+0.05 
"■ '"-0.05 


49.2/27 


7 


0.8lt°- 


, 7C+0.07 


0.29l°°3 


n ifis+0-009 

U.ib»_g ggg 


1 ,0+0.05 
l-l'^-O.OS 


60.5/47 


8 


o.48iE;:i° 


1 71+0.10 

1- '1-0.09 


0.208t°:°?^ 


ior,+0.007 
U.iaU_o gg7 


r, oq+0.04 

U.»3_g g4 


85.5/51 



X-44 



1 


1 10+0. 18 
1-1°-0.20 


2 03+°-°'^ 
^-"'-'-0.04 


39+0-04 
^•■^^-0.03 


164+0-009 

U.lD4_g ggg 


1 26+°-°^ 


35.7/34 


2 


0.83i°:i° 


1 74+0.04 
1- '^-0.04 


030+0.020 

U.OOU_g gig 


190+0-007 

U.iMU_g gg7 


1 90+0.04 
1-^S_0.04 


98.6/87 


3 


0.5l°i 


1 f;o+0.09 


n 3^5+0.05 


,-, 909 + 0.020 
U.ZOZ_g gj^g 


1 40+0.09 

l-4S_g gg 


20.7/22 


4 


0.86l°:« 


-, 01 +0.06 
^-°J^-0.06 


36+0-03 
"■•^"-0.03 


0.198l°:°« 


1 oc;+0.06 
l-'J"'-0.06 


37.6/42 


5 


1 1+0.3 

1-1-0.6 


9 9+0.3 
^•^-0.3 


30+0-06 

U.dU_g g7 


0.11ll°-°f7 


r, 07+0.14 

U.»l'_g.g7 


14.8/15 


6 


75+°-18 

U. '0_g 21 


1 fio+0.09 


34+0-04 
^-■5^-0.03 


214+°-°13 

U.^14_g gi3 


1 40+°-°^^ 
1-^^-0.07 


41.4/34 


7 


0.64l°i« 


1 90+0.08 
l-y'-'-0.07 


99+0.03 
1-^^-0.02 


n ,19+0.008 

U.ll^_g.gg7 


76+0-04 
"■ 'b-0.04 


21.7/29 


8 


0.86l°:i^ 


1 96+°-°° 

l-aD_g g7 


30+0-03 
"■•i"-0.02 


,.,+0.008 

U.141_g ggg 


1 01 +0.04 
l-'Jl-0.04 


37.1/46 



Note. — Col. (1): Observation index corresponding to Column (1) of Tablc|8] Col. (2): Absorption column density in units of 102i cm~2 
Col. (3): Photon index of the powerlaw model. Col. (4): Normalization of the powerlaw model at 1 keV in units of 10~* photons cm^2 s^i. 
Col. (5): Absorbed flux in 0.3—10 keV in units of 10~^2 ergs~i cm~2 s~i. Col. (6): Unabsorbed luminosity in 0.3—10 keV in units of 
10*0 ergs s~^. Col. (7): Best-fit a-nd degrees of freedom. All errors are at 1 ct level. 



(|Kubota fc Makishimal[200l . the maximum i?in in such 
a phase corresponds to the innermost circular stable or- 
bit (ISCO) of the accretion disk, -which is 9 km for a 
non-spinning black hole of one solar mass. This places a 
lo-wer limit on the compact obje ct mass of M > 11 Mq. 
Adopting Equation (9) in Maki shima et al.l (|200Cl[ ). sub- 
stituting ^ = 0.41 and Teg- = T^^/k ("where ^ is the radius 
correction for the maximum temperature, and k is the 
hardening correction -which is equal to / in our paper), 
we have L = 6.0 x 10^'^ (M / Mq)^ (T^s /keV)^ ergs s'^ for 
a non-spinning black hole. Fitting to the data of NGC 
5204 X-1 gives a best-fit relation Lx = (7.9 ± 0.8) x 



10'*^(rcff/keV)'' ergs s^^, and, thus, a black hole mass of 
360 ± 120 Mq. The luminosity, disk size and tempera- 
ture all point to a consistent result that NGC 5204 X-1 
harbors an IMBH -with a fe-w hundred solar masses. We 
note that the derived mass will be larger for a Kerr black 
hole, which has a smaller ISCO and hotter disk than a 
non-spinning black hole. We also note that there is a 
suggestion of high energy spectral curvature in the spec- 
trum indicative of an optically thick corona, which could 
hide the innermost accretion disk leading to an overesti- 
mate of the black hole mass. However, the evidence for 
curvature is weak in the XMM band, and is only seen in a 
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Lx versus F of Antennae X-11, X-16, X-42, and X-44 derived from an absorbed power-law model (data points adopted from 



Fig. 
Table O. 

few bins at energies close to 10 keV. Future observations 
in the hard X-ray band could address this question. 

In Holmberg II X-1, the four data points are not in- 
consistent with L <x. T'^ relation. However, more obser- 
vations are required to accurately measure the spectral 
evolution and test the L oc T** relation before concrete 
conclusions can be made. 

Interstellar absorption and the cool disk model are 
competing in the same energy range. As demonstrated in 
Figure [21 lower abundances will not change the correla- 
tion found for the disk model, but will mildly flatten the 
luminosity and temperature evolution pattern. Higher 
abundances will affect in the opposite way. Therefore, 
for systems like NGC 1313 X-2 where an anti-correlation 
was found between the luminosity and temperature of 
the cool disk, changes of the abundance will not help re- 
cover the L oc T"* relation. For systems like NGC 5204 
X-1 in which the correlation is found, the abundance is 
important in determining the power-law slope between L 
and T. For NGC 5204 X-1, near-sola r abundance was de - 
tcrmincd using XMAf observations ([Winter et a.l 120071 ). 
For Holmberg II X-1, differe nt abundances ar e measured 
using X MMRGS (0.6 solar : [Goad et al.l [2006l) and EPIC 
(1 solar; [Winter et a.ll2007l ) data. If the 0.6 solar abun- 
dance is true, the slope between L and T will be smaller, 
which makes it more similar to NGC 5204 X-1. 

The cool disk interpretation of the spectrum of IC 342 



X-1 seems to be incorrect because the luminosity ver- 
sus temperature data are strongly inconsistent with the 
L oc relation (Figure [7]). For the two observations 
with low luminosities, the confidence levels found by F- 
test of the need for a cool disk component in addition 
to the power-law component are 0.77 and 0.9992, respec- 
tively for observation 1 and 3. Therefore, the cool disk 
component in the 3rd observation is significant at 3.3cr, 
which firms the anti-correlation between the luminosity 
and temperature, but more observations are needed to 
test the results. This low/hot to high/cool evolution de- 
rived from the cool disk model is quite similar to the 
behavior of NGC 1313 X-2 ([Feng fc Kaaret[[2007a[) . The 
best-fit power-law exponent between Lx and Tin is also 
similar for the two systems: -2. 8 ± 0.7 for IC 342 X-1 
and -3.7±0.7 for NGC 1313 X-2. [Poutanen et"all (pOOTi ) 
predicted an anti-correlation between Lx and Tin in a 
model that assumes ULXs are super-Eddington sources 
similar to SS 433 and the soft excess is due to thermal 
emission at the spherization radius where outflows start. 
Both IC 342 X-1 and NGC 1313 X-2 are surrounded by 
super nova remnant like nebulae ([PakuU fc Mirionl[2002l . 
[200l . which may be produced by strong outflows from 
the central source, like W 50 around SS 433. Therefore, 
IC 342 X-1 and NGC 1313 X-2 may represent a class of 
sources with super-Eddington accretion. 



Spectral Evolution of ULXs 



15 



3.2. The Power-law/Comptonization Component 

A power-law model usually provides good fits to the 
high energy part of XMM spectra of ULX, leading to 
interpretation of the emission as from an optically thin 
corona. However, in some cases, spectral curvature is 
evident and a Comptonization mo del, or equ ivalently a 
cutoff power-law mo del, is favored ()Feng fc Ka arct 200f3; 
IStobbart et al.l 12006 ). Fitting with the Comptonization 
model requires an optically thick corona in some or all 
observations of NGC 5204 X-1, Holmberg II X-1, and 
IC 342 X-1. The detected spectral curvature starts at 
energies close to the upper bound of the XMM effective 
energy range. Therefore, non-detection of an optically 
thick corona could be a result of the limited energy range, 
and even with a detection, the significance is usually low 
due to lack of information above 10 keV. The dominant 
luminosity of the power-law or Comptonization compo- 
nent in the XMM band suggests that coronal emission 
is dominant in ULXs. The coronal temperature in some 
ULXs is much lower than in Galactic black hole binaries, 
which could be caused by a cooler and more luminous 
accretion disk in ULX s that seeds and cool s the corona 
(|Done fc Kubotal[200l[Stobbart et al][2006l) . 

3.3. The Thermal Dominant State 

In Galactic black hole candidates, the X-ray emis- 
sion in the thermal dominant state is dominated by a 
geometrically thin optically thick accretion disk, which 
can be well described by the st andard accretion disk 
model (S hakura fc SunvaevI [l973h . The L (x T'^ rela- 
tion in this state is we ll defined by observations (cf. 
iGierlihski fc Don3 l2004f ) and is, perhaps, the best un- 
derstood emission state in black hole binaries. How- 
ever, the thermal dom inant state rarely appears in ULXs. 
iKubota et"aI1 ()2001 ^ found that the ASCA spectrum of 
IC 342 X-2 in the high state could be adequately fitted 
by a standard accretion disk model, implying that the 
source was in the thermal dominant state. This is con- 
firmed by the XMM data in a similar energy range of 
1-10 keV. However, as shown in Figure [3 no correlation 
is found between the luminosity and temperature. The 
four data points appear at two luminosity levels varying 
by a factor of 4. Fitting with an L (xT'^ relation results 
in = 120.9 with 3 degrees of freedom, strongly ruling 
out the model. 

3.4. The Hard State 

The four ULXs in the Antennae and IC 342 X-1 dis- 
played a hard power-law photon index in all or some of 
the XMM observations. Antennae X-16 may be the most 
interesting source among them due to its extremely hard 
spectrum with F = 1.1 — 1.4, which are consiste nt with 
previous Chandra results (|Feng fc Kaaretl[2006b[ ). Some 
Galactic objects with such a hard spectrum are usually 
highly absorbed hard X-ray objects with wind accretion 
from a supergiant companion. However, X-16 has a small 
absorption column density beyond the Milky Way; A'h 
is pegged to the Galactic value in four out of the eight 
observations. The source flux is highly variable, while its 
spectrum is constantly hard and the absorption is low. 
Lack of radio emission rules out that the source emis- 
sion is due to re lativistic bea ming along the line of sight 
such as blazars (jZezas et al.ii2002b) . An analogue to X- 
16 is CXO J024238.9-000055 in NGC 1068, which also 



shows an extremely har d spectrum with little absorption 
(jSmith fc Wilsonl[2003h . 

Antennae X-11 and X-42 have a constant F around 
1.6-1.8, except for one data point for X-42 in obser- 
vation 4 with F — 2.21^q'5i that is softer than most 
others. Antennae X-44 displays a little trend of soft- 
ening at low luminosities until F = 2.0. IC 342 X- 
1 shows a flux change by a factor of a few with F < 
2.1. These behaviors are consistent with the classifica- 
tion as in the hard state with F = 1.5 — 2 .1 in black 
hole binaries (|McClintock fc Remiflardl I2006D . The two 
brightest (in observed flux) ULXs in nearby galaxies 
are X 41. 4-1-60 and X42. 3-1-5 9 in the starburst galaxy 
M82 (jFeng fc Kaaretl l2007bl : iKaaret et all 120081) . and 
also show hard power-law spectra. X41. 4-1-60 remains 
in the hard s tate even at the highest fluxes observed 
(jKaaret et al.l 12008). The Antennae galaxies and M82 
both have high star formation ra tes, and the ULXs in 
them share similar properties. iBerghea et al.l (j2008l ) 
found ULXs tends to be hard when the luminosity in- 
creases. The hard state represents a luminosity of a few 
or as high as 30 percent of the Eddington hmit (e.g., 
GX 339-4: iMivakawa et al.l[20M l. Therefore, the mass 
of the compact object in these ULXs, if they are in the 
hard state, is at least a few hundred solar masses. 

3.5. The L — F Correlated Phase And Steep Power-law 

State 

NGC 5204 X-1 was probably in an intermediate state 
with hybrid properties from the thermal dominant and 
steep power-law states during those XMM observations. 
In a 50 ks Chandra observation made in 2003, the source 
made a transition to the steep power-law state with F = 
2.6 or 2.9 (dep ending on models) a t a luminosity of a few 
Iq39 gj-gg g-i f Rob erts et al.ll2006h . This certainly broke 
the tight correlation found between Lx and F as shown 
in Figure [1] and is very similar to the behavior of NGC 
1313 X-1 iFeng fc Kaaretl [2006ah . NGC 1313 X-1 also 
showed a tight correlation between Lx and F at F < 2.4, 
and then jumped into the steep power-law state with 
F = 3.1. Holmberg II X-1 might also have undergone 
a similar transition from an Lx — F correlated phase to 
a steep power-law state (see Figure [5|) based on the four 
observations. Another similarity between Holmberg II X- 
1 and NGC 5204 X-1 is that they both show correlation 
between Lx and Tin. This analogy suggests that soft 
excesses in NGC 1313 X-1 may be due to disk emission 
around an IMBH and may follow an L cx evolution. 

This L — F correlated phase could be common in ULXs 
with an IMBH accretor. It shows hybrid properties from 
the thermal dominant state (with an i oc T** relation), 
and the steep power-law state (with strong Comptoniza- 
tion). Sources in this phase show strong variability and 
occasionally make a transition to the steep power-law 
state. 

3.6. IC 342 X-2 And The Peculiar Soft Feature 

The source shows a peculiar feature at energies below 
1 keV. A collisionally-ionized diffuse gas model (apec) 
could adequately fit this feature. However, we argue 
that the best-fit parameters derived from the apec model 
are unphysical. The normalization of the apec model is 
Na = IQ^^^nlV/ {4ttD'^), where is the electron density 
in cm^'^, V is the volume in cm'^, and D is the distance to 
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the source in cm. The total kinetic energy of the plasma 
is = ikTn^V = 6 X lO^kTNAD^/ue ergs s'^, where 
kT is the plasma temperature in keV. Substituting kT, 
Na, and D, we obtain E ^ lO^^/rie or lO^^/ue ergs s~^ 
for the two observations. These energies are too high for 
an accreting system even with an IMBH accretor in an 
extremely dense environment. The variation of the apec 
component also excludes the possibility that the soft ex- 
cesses could be caused by a supernova remnant or other 
large scale emitter. We tried other models but none was 
able to fit the soft excesses as well as a thermal plasma 
model. The nature of soft spectral feature in IC 342 X-2 
is puzzling. 

3.7. Timing Properties 

The lightcurve of Holmberg II X-1 showed variability 
in the energy range of 1-4.5 keV (see Figure (5]). Pe- 
riodograms were calculated and a major peak around 
20 ks was remarkable, which corresponds to multiple 
spikes in the lightcurve at the same timescale. It is 
unclear whether or not the variation is periodic; more 
observations are needed to test that. If this is a peri- 
odic variation and modulated by the orbital motion of 
the binary, then the observed timescale could be the 
orbital period. The density of the companion star can 
thus be inferred under this scenario assuming Roche- 
lobe overflow as p = llOPj^.^ = 3.6 (g cm^'^), corre- 
sponding to a low mass late type main sequence star 
(e.g., MO V). Unevolved low mass companions have not 
yet been discovered in black hole binaries, but exist in 
ultra-compact binaries with a neutron star accretor. It 
is inconsistent with the general picture that ULXs in 
star forming regions are usually high mass X-ray bina- 
ries (e.g., Feng & Kaaret 2006b). However, ULXs in el- 
liptical galaxies, especially those associated with globu- 
lar clusters, must be low mass X-ray binaries. On the 
other hand, if the variability is quasi-periodic oscilla- 
tions, it has a central frequency much lower than those 
found in other ULXs (Strohma ver fc Mushotzkvl 120031 : 
IStrohmaver et all 12007) and would be a good argument 
for an IMBH. Therefore, interesting consequences are ex- 
pected if the nature of this 20 ks variability are revealed 
with future observations. 

We report no other significant short term timing noise 
from these observations due to the limited sensitivity, 
though it is expected in the hard state and sometimes in 
the steep-power law state. 

4. CONCLUSION 



In this paper, we analyzed the spectral evolution of 
NGC 5204 X-1, Holmberg II X-1, IC 342 X-1 and X-2, 
the Antennae X-11, X-16, X-42, and X-44. New results 
obtained in this paper and those reported in the litera- 
ture (including ULXs in M82 and NGC 1313) were dis- 
cussed together. The main conclusion of the paper is to 
report two major spectral classes in ULXs: one is the 
hard state and the other is the L — T correlated phase. 

Many ULXs, including sources in M82 and the An- 
tennae, are found in the hard state, in which the source 
shows strong flux variability but the spectral index re- 
mains hard and constant. From analogy with the be- 
havior of Galactic black hole binaries in the hard state, 
the emission of these ULXs should be sub-Eddington and 
their masses are estimated to be at least a few hundred 
solar masses. The Antennae X-16 is an extremely hard 
ULX, with a power-law photon index F = 1.1 — 1.4; 
sources with such a hard spectrum, low absorption, and 
high luminosity are rare and difficult to understand. 

The second class of ULXs show behaviors including: 
(i) the luminosity is positively correlated with the power- 
law photon index, (ii) the source sometimes breaks the 
correlation and jumps into the steep power-law state, 
and (iii) the luminosity seems to vary with the 4th order 
of the disk temperature indicating that the emission is 
from an accretion disk. Sources in this class includes 
NGC 5204 X-1, Holmberg II X-1, and NGC 1313 X-1. 
We interpret these behaviors as a hybrid between the 
thermal-dominant and steep power-law states. Sources 
in this class are thought to contain compact objects more 
massive than stellar-mass black holes inferred from their 
luminosity versus disk temperature pattern. 

Some sources, like IC 342 X-1 and NGC 1313 X-2, 
show a negative correlation between the luminosity and 
disk temperature. If the negative correlation can be con- 
firmed with future observations and commonly found in 
other ULXs, they may represent a third class of ULXs 
with super-Eddington emission. 

Also, we report a peculiar soft excess in IC 342 X-2, 
which is hard to explain, and strong variation on the 
timescale of 20 ks from Holmberg II X-1. More observa- 
tions are needed to test these results. 
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